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Abstract 



, We calculate the production cross sections for D*~^D*~ , D~^D*~ and D~^D~ in e~^e~ 

' annihilation through one virtual photon in the framework of perturbative QCD with 

constituent quarks. The calculated cross sections for D*~^D*~ and D^D*~ production 
. are roughly in agreement with the recent Belle data. The helicity decomposition for 

I D* meson production is also calculated. The fraction of the D'^D^ final state in 

^f) • e'^e~ — > D*~^D*~ process is found to be 65%. The fraction of DD^ production is 

100% and DD^ is forbidden in e~^e~ annihilation through one virtual photon. We 
further consider e~^e~ annihilation through two virtual photons, and then find the 
O ' fraction of DD^ in e~'"e~ DD* process to be about 91%. 

^ i PACS number(s): 12.40.Nn, 13.85.Ni, 14.40.Gx 

p • Recently, heavy meson production in e^e~ annihilation has become a very interesting 

^ ■ subject both experimentally and theoretically. For instance, for charmonium production in 
e^e~ annihilation at the B-factory energy s/s = lO.GGeV, there are large differences be- 
tween the experimental data [TJ |2l E] and the calculated cross sections for both exclusive 
^ I processes [3] and inclusive processes jS]. Even by including the effects of e~^e~ annihilation 
' into two photons in the exclusive double-charmonium production [Hj and inclusive char- 
monium production [7j, the large discrepancies still exist [S]. Moreover, most recently the 
Belle Collaboration has measured jH] the charmed meson pair D^*"^^ D^*^~ production in 
e~^e~ annihilation and also found a large differences between observed data and theoretical 
predictions PU]. The measured cross sections are a{e^e~ D*'^D*~) = 0.65 ±0.04 ±0.07 pb 
and a{e~^e~ — D^D*") = 0.71 ± 0.05 ± 0.09 pb, which are, however, lower than those 
predicted in Ref.jinj by an order of magnitude (see [Hj for the comparison). The Belle mea- 
surements of the exclusive D meson pair production in e~^e~ annihilation could be another 
challenge to the theoretical studies. In order to understand these production processes, in this 
letter we will present a calculation in the framework of perturbative QCD with constituent 
quarks. Namely, we will treat the charmed mesons D^*^ as bound states of a constituent 
charmed quark and a constituent light antiquark, and the virtual photon will couple to the 
charm quark, which will subsequently emit a light quark pair with constituent quark masses 
through a virtual gluon (see the upper diagram in Fig. Q). The virtuality of the virtual gluon 
could be large enough for the application of perturbative QCD. In addition, the case of the 
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virtual photon coupled to the light quark will also be considered (see the lower diagram in 
Fig.^. In the following we will report our calculated results for both the total cross sections 
and the helicity decomposed cross sections for these charmed meson pair production in e^e~ 
annihilation at y/s = 10.6GeV, . 

Following the method in Ref. [TTj, the amplitude for producing the heavy quark (Q) 
-light antiquark (g) bound-state (Qq) is given by 

MP)= E f -^,^L,LA^){LLz;SSz\JJz)M{P,k), (1) 

Lz,Sz-^ ^^'^> 

where 

M{P,k) = OrTssz{P.k), (2) 

and Oy represents the short-distance interaction producing the Q and g in a specific bound- 
state, k is the relative momentum between Q and q. \l/i^L^(k) is the Bethe-Salpeter bound 
state wave-function. Vsszi.Pi up to second order in k, is given by 



TssziP. k) = + 1^ I SSzMrP - k, s)u{fP + fc, s), 

^V'^'^ + ^l 2m, )\-^{P.Sz))\ 2mg ) 

where 7^ is for spin-singlet state and — /(P, 5^) is for spin-triplet state. 

For the S-wave state, we have M{P, k) = M{P, 0). After integrating over k, the amplitude 
A{P) is related to the origin of the radial wave-function Rs{0) 

'''^*„o^^. (4) 



Using Eq. (PIH), the calculation of the cross section for the process showed in Fig. [T]is 
straightforward. We get the cross section as follows 



da= ' ^ ' \ l \M\^dx, (5) 

8s V s 

where x = cos 6', 6 is the angle between pi and ps. The meson mass m equals rric + md in 
the leading oeder non-relativistic approximation. For the e~^e^ D*~^D*~ process, | M p 
reads 



128(mc + ■mdy{4ml + Smcma + 4:ml — s) 

X [-48mf - 192ml^md - 288mJ°m^ - 192m^m^ - 96mlmj - \92m\m\ 
-288m|;m^ - 192mlml - QOmjml - 48m^m^ - 72m^m^° - 48mcm^^ 
— 12m^ — Sml^s — ASmlmdS — 80m^m^s — 64m^m^s — 56m^m^s 
— 64m^m^s — bOm^m^s — 28m^m^s — 20m^m^s — 12mcm^s — 2m}^s 
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+192m^mj] + 96m^m^ + 192m^m^ + 288m^m^ + 192m^m^ + 60m^m^ 
+48m^m)^ + 72m^m^° + ASmcTJiY + I2mf - 8mJ°s - 16m^mdS - 16m^m^s 

a A f^f^ Ad QT 

— 16m^m^s — 24m°m^s — 32m^mj^s — 18m^m^s — 4mj:m^s 

—Amlm^s — Arricm^s — 2ml^s + 4m^m^s^ + Am^m'^s^ + m^m^s^)x^]. (6) 



For the e+e ^ D^D* process, one should replace y 1 ^ with 

in the phase space integration. | M p is given by 

For the e~^e^ D^D^ process, | M p is as follows 



^ 2 128(mc + mdY{Aml + Snicnid + 4m^ — — x^) 

X (4m^ + Smlnid + Am^ml + 2m^m^ + 4mcm^ + 2m^ — 2'm\mdS — mcm\s). (8) 
The Mandelstam variables are defined as 



t = {Vz-Pif = m^--\l-\l X , (9) 




M = (P3-P2) =m^-- 1 + W1 + x\. (10) 



We notify that for e^e~ D^D*~ process, the formula of Mandelstam variables t and 
u are a little different from Eq. Q and Eq. (jlUp because of the mass difference between 
and D*~. As the difference is small, the change for the cross section is negligible. For 
simplicity of the module of amplitude, we still use the definition in Eq. (0) and Eq. (fTUI) . 

In the numerical calculation, the input parameters follows 

a = 1/137, a, = 0.26. (11) 

For the Coulomb-plus-linear potential case, the value of the wave function at the origin 
for charmonium can be found e.g. in |12|. With rric = 1.84 GeV, for the S-wave charmonium 
it could be given by 

|i?5(0)P = 1.454 GeVl (12) 

We then use the potential scaling rules to get a rough estimate for wave function at the 
origin for the charmed meson. We fix the constituent c-quark mass at 1.6 GeV. The cross 
sections for !)(*)+£)(*)-, and D+D- are found to be 

a(e+e" ^ 7* ^ D*+ D*-) = 0.532 pb, (13) 
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a{e+e- 7* ^ D+D*-) = 0.699 pb, (14) 

a(e+e- 7* ^ D+D-) = 0.098 pb. (15) 

We also calculate the polarized and D^D^*^~ production. The polarized cross 

section can be calculated by defining the longitudinal polarization vector as follows 

P =J7 , (16 

M n ■ p 

where = and = (1, —p/ \ p\). 

The polarized modules of amplitude for e+e^ — > 7* — ^ £){*)+ jji*)- jg listed in Eq. El 



— 128(mc + mdY{2ml + mj^)^(4m^ + Smcirid + 4m^ — s)(l + 



M \'= ^^^^ ^ ' , \ ' ^ " 2 — (17) 



The corresponding cross sections are 



a(e+e- 7* D^*^^dP^) = 0.347 pb, 



Using the results in Eq. (fT^ and (fTHj) . we get 

cr(e+e- ^ D^*^^dP^) 
a{e+e- £)(*)+£)(*)-) 



65.2% (19) 



For D~^D^*^~ case, the longitudinal cross section for one photon process must be zero, 
since it is forbidden by parity conservation and angular momentum conservation. The effec- 
tive ^DD* vertex in Fig ^ only can have the form as 

~ e''''^''F,,{q)GpM<PiP3), (20) 

where F^y{q) represents the field of the virtual photon, Gpaij)^) represents the field of the 
jji*)^ meson and ^(ps) represents the field of the meson. In the momentum space, the 
effective vertex is as follow 

C^nt-^e^'''"^q,tl{q)p,,4{p,). (21) 
If we choose the center-of-mass frame of the e"'"e~, the formula in Eq. ()21|) becomes 

C^nt ~ e''P"qotl{q)p,p4{p,) 

- go(p4Xe'-''(p4)) ■€-*(?). (22) 

Because the space component of the longitudinal polarization vector of the D* meson is 
parallel to the momentum of the D* meson, from Eq. ()22|) . one knows that the longitudinally 
polarized D* meson cannot contribute to the process listed in Fig. ^ The longitudinal D* 
only comes from two photons process, as listed in Fig. |21 



4 



The modules of amplitude for e+e 27* D~^D^*^ and e+e 27* D^d!^^ are 
listed, respectively, in Eq. (j23|) and Eq. (j24|l . 

I ^ I 8(mc — m^Yijnc + mdYirnl + Snicrrid + m^)^(4(mc + mrf)^(l — x^) + s(l — x^)) 

(23) 

I ^ |2_ 32(mc - nidYinic + mrf)^°(mc2 + "jnicmd + + x^) 

81m^m^s^(2mc + md)^(mc + 2md)'^ 
The corresponding cross sections are 

a(e+e" ^ 27* ^ D+D^*^-) = 0.091 pb (25) 

cT(e+e- ^ 27* ^ D+dP~) = 0.018 pb. (26) 
Then one can get the ratio of cross section for to total cross section as follows 



90.8%. (27) 



In summary, motivated by the measurement of the e"'"e^ [)(*)+[){*)- cross-sections [H], 
we calculate the cross-sections for /)(*)+_D(*)^j D^D^*^" and D^D^ in e~^e~ annihilation 
through one photon. The calculated cross-sections are roughly in agreement with the ex- 
perimental data. The fraction of the D^^-''^ d!^''~ final state in the e+e" /)(*)+£)(*)- 
reaction is also calculated, the ratio of 65% is in some deviation from the Belle data, which 
is (97 ±5)% p|. Moreover, for the D^D'^*'>^ production case, in Ref. 0, it is claimed that 
e+e- ^ D+D^*^- is saturated by Dl final state (95.8 ± 5.6)%, in good agreement with the 
predictions of Ref. fUl- Under our analysis, the final state is forbidden by parity con- 
servation and angular momentum conservation. If the e+e~ 27* D~^D^*^~ process is 
further considered, the fraction of the final state in e~^e~ — >■ D^D^*^^ process is 90.8%. 
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Figure 1: Feynman diagrams for e+e" ^ 7* ^ 
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Figure 2: Feynman diagrams for e+e" 27* £)(*)+£)(*) 
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